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ULTRASONIC  CHARACTERIZATION  OF  HIGHLY 
ATTENUATTVE  THICK  COMPOSITES 


1.  INTRODUCTION 

In  recent  years  there  has  been  tremendous  progress  in  the  development  of  new  kinds  of  structural 
materials.  It  is  now  possible  to  combine  together  several  desirable  material  properties  which  were 
often  mutually  exclusive  to  manufacture  custom  materials  according  to  specific  requirements.  As 
an  example,  fiber-reinforced  composites  can  be  utilized  in  structural  components  which  require 
high  strength  and  light  weight.  Another  reason  for  their  popularity  is  the  possibility  of  customizing 
their  composition  according  to  specific  requirements.  In  fact  one  can  use  different  kinds  of  fibers, 
or  a  different  fiber  layout  and  density,  or  a  different  material  as  a  “matrix”  according  to  the 
intended  use  and  required  properties.  Possibilities  range  from  highly  ordered  patterns  (e.g., 
fiberglass  woven  cloth  layers  imbedded  in  epoxy)  to  randomly  dispersed  chopped  fiber  strands. 

The  variety  of  new  materials  already  developed  and  the  continuing  trend  towards  highly 
specialized  composites  suggest  that  attention  be  devoted  to  the  problem  of  their  characterization. 
In  fact,  safety  and  economic  considerations  require  reliable  techniques  for  NDE  of  these  new 
materials.  It  is  also  necessary  to  gain  a  good  understanding  of  the  effect  of  the  different  constituents 
and  the  role  of  the  interfaces  in  order  to  predict  properties  with  the  goal  of  optimization,  in  view  of 
specific  requirements.  A  severe  difficulty  in  the  NDE  of  many  of  these  newly  developed  materials 
is  their  high  absorptivity  with  respect  to  the  interrogating  energy,  e.g.  ultrasonic  waves.  In  fact, 
from  NDE  point  of  view,  a  material  specimen  can  be  defined  thick  or  thin,  not  according  to  its 
physical  thickness  (as  expressed  e.g.  in  mm)  but  to  its  attenuation  coefficient  for  the  interrogating 
energy.  This  concept  will  be  clarified  further  in  Section  2. 

There  are  a  variety  of  problems  which  arise  in  the  NDE  of  thick  composites,  depending  on  the 
physical  properties  of  the  material  to  be  characterized,  the  nature  of  the  flaws  to  be  investigated 
and  the  special  requirements  for  NDE  specifications.  In  fact  some  NDE  techniques  are  more 
suitable  for  certain  materials  than  others.  The  possible  flaws  that  occur  in  these  materials  can  be 
overwhelming:  there  may  be  porosities,  contaminants,  inclusions,  delaminations,  moisture,  fiber 
misalignment  or  breakage,  matrix  aging  or  inhomogeneities,  etc.  We  wish  to  remark  here  that  the 
term  “flaw”  may  be  misleading  because  some  of  the  “flaws”  may  really  be  desirable  features  such 
as  irregularities  in  the  heterogeneous  composition  of  the  material.  The  signal  due  to  the  above 
detrimental  flaws  in  thick  composites  is  generally  a  small  perturbation  on  the  received  signal  which 
is  in  the  noise  level.  For  the  purposes  of  detection  of  competing  signals  due  to  these  flaws,  one 
requires  well-characterized  high  power  transmitters  and  sensitive,  high  resolution  receivers  with 
appropriate  signal  processing.  Among  the  “special  requirements”  which  might  affect  the  selection 
of  a  particular  NDE  technique  over  others,  we  may  mention  the  availability  and  cost  of  the 
apparatus,  its  ease  of  use  in  the  field,  the  time  required  for  the  job  and  the  accessibility  of  the 
component  to  be  inspected,  as  the  deciding  factor.  A  consideration  which  may  prove  to  be  of 
primary  importance  is  whether  the  detection,  location  and  characterization  of  the  flaws  in  the 
component  needs  to  be  performed  during  the  processing,  fabrication  or  service  phase  of  the  work. 
Another  important  consideration  is  whether  the  data  presentation  and/or  accept-reject  criterion  is 
qualitative  or  quantitative.  In  many  cases  the  detection  and/or  evaluation  of  the  severity  of  the  flaw 
may  be  greatly  eased  by  a  complete  automation  of  the  NDE  process. 
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It  is  important  to  assess  and  predict  nondestructive^  the  presence  of  irregularities,  porosities, 
inclusions,  and  flaws,  such  as  cracks  and  delaminations.  For  example,  for  composite  materials  such 
as  fiber-reinforced  polymers,  it  may  be  necessary  to  map  ultrasonically  the  orientation  and 
alignment  of  the  fibers  and  to  detect  irregularities  such  as  misalignment  or  heterogeneities  of  the 
fibers.  In  each  case,  one  needs  to  measure  variations  in  ultrasonic  attenuation  and/or  phase  velocity 
to  assess  the  state  of  the  material. 

A  large  number  of  papers[l-19],  applicable  to  the  problem  of  characterizing  composite  materials 
already  exist.  In  this  report  (Section  3),  we  shall  limit  ourselves  to  various  ultrasonic  techniques 
which  we  developed  for  application  to  the  characterization  of  highly  attenuative  thick  composites. 


2.  NDE  DEFINITION  OF  THICK  COMPOSITES 

Specimens  made  of  composite  materials  are  usually  thin  in  one  or  more  dimensions  (e.g.  a  few  mm 
or  a  few  cm  thick).  If  the  specimen  is  not  of  uniform  thickness,  the  term  “thick”  generally  refers  to 
the  load-bearing  thick  section.  However  from  an  NDE  point  of  view  what  matters  is  the  attenuating 
power  of  the  specimen,  with  respect  to  the  particular  kind  of  interrogating  energy  used.  Therefore 
the  thickness  of  the  material  which  can  be  inspected  by  NDE  is  restricted  by  the  signal-to-noise 
ratio.  This  can  be  better  understood  with  a  few  examples. 

The  amplitude  of  an  acoustic  wave  propagating  in  a  medium  in  the  z -direction  can  be  written  as 

A(z)  =  A0  e*1  F(kz-<ot)  (2.1) 

where  the  attenuation  coefficient  a  (in  nepers/cm  or  nepers/inch)  depends  on  the  material 
properties.  A0  is  amplitude  at  z  =  0.  F  is  a  function  of  wave  vector  k,  and  angular  wave  frequency 
to  and  defines  the  shape  of  the  interrogating  ultrasonic  energy. 

Assuming  that  the  noise  level  amplitude  is  given  by  AHl,  the  equation 

4  (*„„)  =  A„,  (2.2) 

defines  the  maximum  distance  zmax  for  which  a  determination  of  A  (z)  can  be  relevant  i.e.,  for  the 
ultrasonic  NDE  of  a  material,  the  thickness  d  must  be  such  that,  for  a  given  frequency,  the 
amplitude  of  the  transmitted  signal  A  (d)  >  AHl.  If  A(d)  <  Aml  the  material  is  defined  as  thick. 
Therefore,  the  decision  of  whether  a  material  is  thick  or  thin  for  the  purposes  of  ultrasonic  NDE 
depends  on  the  acoustic  properties  of  the  material.  For  some  highly  heterogeneous  composites 

a  =  14.39  nepers/inch  =  125  dB/inch  (2.3) 

which  means  that,  if  input  signal  A0  =  1  v,  the  amplitude  of  the  electric  signal  representing  the 
acoustic  wave  on  the  oscilloscope,  after  a  penetration  of  one  inch  is  given  by 
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A(l")=0.56  F(kz-<at)  pV 


(2.4) 


which  could  be  below  the  noise-level  of  many  ultrasonic  systems.  However,  if  instead  of 
ultrasound,  a  different  type  of  interrogating  energy,  e.g.  electromagnetic  waves  is  used  for  NDE, 
the  criterion  for  inspectable  thickness  would  depend  upon  the  electromagnetic  properties  of  the 
material.  All  types  of  composites  can  be  inspected  by  ultrasound  whereas  electromagnetic  NDE 
techniques  [20-22]  are  applicable  only  to  the  composites  which  are  electrically  conductive  such  as 
graphite/epoxy  and  may  be  able  to  inspect  greater  thickness.  Moreover,  these  composites  also  do 
not  have  high  acoustic  absorption  coefficient  as  compared,  e.g.,  to  kevlar/epoxy  composites  and 
thus,  do  not  present  S/N  ratio  problem.  We  restrict  our  ultrasonic  study  the  latter  class  of 
composites,  some  which,  depending  upon  the  fiber  layout,  can  have  very  high  attenuation 
coefficient. 


3.  ULTRASONIC  METHODS* 


Ultrasonic  techniques  offer  several  features  which  in  many  cases  make  them  particularly 
convenient  for  NDE.  They  arc  safe,  cost-effective  and  easy  to  use  in  different  situations  like  in  the 
laboratory  or  in  the  field.  The  high  ultrasonic  reflectivity  of  narrow  defects  such  as  cracks  and  their 
accessibility  to  ultrasound,  even  deep  into  the  material  may  be  an  important  factor  for  choosing  an 
ultrasonic  inspection  over  other  NDE  techniques.  Interaction  of  these  waves  with  the  elastic 
properties  of  the  composite  can  be  used  to  detect  and  evaluate  flaws. 

Attenuation  of  ultrasonic  waves  propagating  through  a  material  is  an  important  physical  parameter. 
It  gives  us  an  indirect  measure  of  the  composition  of  the  material.  Certain  constituents,  such  as 
fibers  in  the  host  matrix,  can  be  added  intentionally  in  order  to  create  materials  such  as  composites 
which  are  lighter  and  stronger  than  conventional  structural  materials,  such  as  steel,  aluminum  and 
titanium,  etc.  In  order  to  characterize  these  composite  materials  nondestructively,  one  needs  to 
measure  the  attenuation  and/or  velocity  of  the  propagating  ultrasound  through  such  materials  and 
relate  these  acoustic  parameters  with  the  material  properties  or  variation  thereof.  In  this  paper  we 
concentrate  on  the  measurement  of  attenuation  in  materials  which  are  highly  attenuative. 

Ultrasonic  pulse-echo  systems  have  been  widely  used  in  recent  years  to  measure  the  attenuation 
and  ultrasonic  group  velocity  with  applications  for  tissue  characterization  in  medicine  [23-26]  and 
material  characterization  in  NDE  [27-30].  Pulse-echo  techniques  require  plane  parallel  specimens 
of  low  attenuation  so  that  multiple  echoes  can  be  set  up  in  the  specimen.  They  do  not  measure 
absolute  attenuation  because  of  uncertainties  due  to  reflection,  transmission,  and  surface  scatter¬ 
ing,  although  losses  due  to  these  effects  can  be  partially  taken  into  account  by  processing  multiple 
echoes.  Due  to  a  lack  of  spectral  purity,  they  do  not  provide  us  with  phase  velocity,  which  is  a 
more  relevant  quantity  for  the  determination  of  elastic  constants,  stresses,  texture,  and  material 
phases  in  the  material.  Ultrasonic  attenuation  due  to  impedance  mismatch,  surface  roughness,  and 
diffraction  is  generally  irrelevant  as  far  as  the  material  characterization  is  concerned.  The  more 
relevant  quantities  are  the  attenuation  due  to  absorption  aabs,  and  to  scattering  ajc.  Thus  accurate 
and  absolute  values  of  the  attenuation  alol  =  aaki  +  atc  and  phase  velocity  can  characterize  a  mac- 
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roscopically  homogeneous  material  almost  completely. 

The  difficulty  of  characterizing  thick  components  lies  in  the  fact  that  the  detected  signal  is  often 
masked  by  noise.  It  is  not  possible  to  observe  multiple  echoes  in  a  plane  parallel  specimen  due  to 
high  attenuation.  Consequendy  one  cannot  determine  the  attenuation  coefficient  from  the  slope  of 
the  envelope  of  these  multiple  echoes.  To  overcome  this  difficulty,  we  developed  various 
transmission  techniques,  as  described  in  this  section,  to  characterize  these  composites. 


3.1  Pulse-Transmission  Technique 

There  are  several  methods  available  for  measuring  the  attenuation  [31].  However,  none  of  these 
techniques  is  suitable  for  highly  attenuative  materials.  We  demonstrate  how  the  above  difficulties 
of  poor  S/N  can  be  circumvented  by  use  of  a  ultrasonic  pulse  transmission  technique  in  conjunction 
with  a  suitable  choice  of  specimen  geometry,  such  as  steps  or  wedges  (Fig.  1). 

Whenever  ultrasonic  waves  (elastic  waves  of  high  frequency)  are  propagated  through  a  material, 
they  lose  energy  through  absorption  and  scattering.  Scattering  from  the  defects  redistributes  energy 
away  from  the  direction  of  propagation  while  absorption  converts  the  ultrasound  energy  into  heat. 
Because  of  the  loss  of  energy  (attenuation),  the  amplitude,  A ,  of  a  plane  wave  propagating  in  the 
z-direction  can  be  written  as: 


A  -  A^e*11  expli(2nft-kz)]  (3.1.1) 

where  /  is  the  frequency,  *  =  y  is  the  wave  number,  and  o  is  the  total  attenuation  factor  (usually 

expressed  in  dB/cm).  In  general,  a  is  a  function  of  frequency  and  increases  with  frequency.  Use  of 
low  frequency  may  increase  the  penetration  depth  but  would  limit  us  to  poor  resolution. 

Equation  (3.1.1)  leads  to  an  expression  for  determination  of  a  as 


a  = 


1  (Z|)  ~ 

l*i-*d  "Uj<*2>  - 


nepers/cm 


(3.1.2) 


A 

Recalling  that  A  ( dB )  =  201og— ,we  can  rewrite  a  as 

a  =  - — !_[A.(dB)  -A,  (dB)  ]  dB/cm  (3.1.3) 

l*i -*al 

For  this  work  we  used  specimens  made  of  pure  epoxy  and  Fiber-reinforced  composites.  These 
composites  consisted  of  fiberglass  woven  cloth  layers  imbedded  in  epoxy.  There  are  also  chopped 
fiberglass  strands  dispersed  randomly  in  the  epoxy  matrix.  One  composite  specimen  consisted  of 
a  14x7  cm  plate  with  three  equal  steps  of  thickness,  z,  =  0.64 cm,  z2  -  l.Olcm  and  z3  =  1.53 cm 
respectively.  The  other  specimen  consisted  of  an  8°  right  angle  wedge  with  a  base  14.5  x  7cm.  A 
third  specimen,  in  the  shape  of  a  wedge  of  13  8° ,  was  made  out  of  epoxy.  This  specimen  was  used 
to  check  the  validity  of  the  experimental  set-up  (Fig.  2).  A  broadband  pulse  is  transmitted  normal 
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to  the  (x,  y)  face  of  a  properly  designed  wedge  shape  specimen.  Another  property  aligned 
transducer  receives  the  transmitted  signal.  The  transmitting  and  receiving  transducers  are  fixed 
relative  to  each  other  and  scan  together  in  the  x-direction. 

Fig.  3  shows  the  attenuation  curve  (dashed  curve)  for  the  first  two  steps  (z,  -  0.64cm  and 
z2  =  l.Olcm)  of  the  composite  specimen.  We  assume  that  the  attenuation  due  to  water  is  negligible. 
By  averaging  the  amplitudes  A,  (z,)  and  A2  (z2) ,  one  obtains  A,  -  -IdB  and  A2  -  -12 dB,  which  yield 
a«  31  dB/cm.  The  oscillations  in  the  attenuation  are  due  to  variations  in  the  fiber  concentration. 

Regions  of  lowest  fibers  concentration  (mostly  matrix)  yield  maxima  and  regions  of  the  highest 
fiber  concentration  correspond  to  minima  in  the  amplitude.  By  averaging  the  maxima  and  minima, 
one  obtains  aJO- 19  dB/cm  and  a^-46  dB/cm.  If  the  fibers  density  were,  at  least  locally,  low,  then 
alo  would  approximate  the  attenuation  coefficient  for  the  matrix.  However,  randomly 
interdispersed  chopped  fibers  in  the  matrix  lead  to  deviation  from  this  assumption.  Attenuation 
curves  for  different  values  of  y  can  be  used  to  map  the  fibers  layout  and  detect  the  regularities  (lack 
of)  in  the  fibers  pattern  and  concentration. 

To  illustrate  the  utility  of  wedges  for  measurement  of  attenuation,  we  used  an  epoxy  wedge 
specimen.  The  advantage  of  using  the  wedge  is  that  the  path  of  propagation  can  be  changed 
continuously  (instead  of  finite  steps)  from  z=  0  to  zmtI  at  which  the  received  signal  is  in  the  noise 
level.  Since  z  =  xtanO,  for  a  fixed  6,  the  path  of  propagation  z  can  be  varied  continuously  by 
scanning  the  pair  of  transducers  along  the  x-axis.  The  amplitude  of  the  received  signal  as  a  function 
of  z  for  a  pure  epoxy  specimen  is  shown  in  Fig.  4a.  (The  scale  on  the  left  is  nonlinear  because  the 
data  system  plots  the  data  in  this  fashion).  The  curve  is  very  nearly  exponential,  so  the  attenuation 
can  be  easily  determined  as  its  average  slope  (See  Eq.  3.2.1).  This  yields  a  -  9 dB/cm.  Figure  4b 
shows  the  perturbation  on  the  attenuation  due  to  scattering  from  a  simulated  flaw  in  this  specimen. 
This  perturbation  on  the  attenuation  curve(exponential)  expected  for  an  unflawed  specimen  can  be 
used  as  a  direct  quantitative  measure  of  the  signature  of  the  flaw. 

Figure  3a  shows  how  the  *bove  mentioned  technique  can  be  used  to  study  a  composite  reinforced 
with  tweed-like  woven  fiberglass  lamina.  As  the  transducers  scan  the  specimen,  the  propagating 
ultrasound  is  intercepted  alternatively  by  regions  of  cross-weave  and  the  matrix.  Due  to  strong 
scattering  by  cross-weave  regions,  the  propagating  ultrasouno  is  highly  attenuated  and  the  re¬ 
ceived  signal  drops  sharply.  The  oscillations  show  the  periodicity  of  the  cross-weave  of  the  lamina 
in  the  matrix.  This  has  a  spatial  period  of  approximately  6.5  mm  over  most  of  the  specimen.  This 
agrees  well  with  the  periodicity  of  the  fiber  weave  measured  optically.  The  exponential  envelope 
(due  to  change  in  path  of  propagation)  give  the  attenuation,  which,  in  the  region  of  cross-weave, 
is  -50  dB/cm  (lower  envelope)  and  in  the  region  of  the  matrix  is  -30  dB/cm  (upper  envelope).  If 
there  is  any  misalignment  or  in-  homogeneity  or  gross  local  flaws,  the  attenuation  curve  will  have 
local  deviation  from  the  exponential  and,  in  the  case  of  the  composite  specimen,  lack  of  periodic¬ 
ity  as  is  evident  from  Fig.  5b.  For  example  the  region  z  =  0.3  to  0.5cm  (at  x  =  1.8  to  3.6cm)  is 
flawed  for  a  particular  value  of  y,  when  compared  with  a  similar  region  for  a  different  value  of  y 
(Fig.  5a).  Thus  this  technique  can  be  used  to  map  the  corresponding  irregularities  of  the  fibers  and 
matrix  flaws. 

Thus  we  have  demonstrated  that  it  is  possible  to  measure  directly  the  attenuation  in  highly  absorb- 
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ing  materials  by  using  a  transmission  technique  and  a  properly  designed  specimen.  The  internal 
composition  of  the  material  is  clearly  depicted  by  the  attenuation  curve. 

32  Ultrasonic  Interferometric  Technique 

In  this  section  we  will  discuss  tow  one  can  measure  both  attenuation  and  phase  velocity  simulta¬ 
neously  in  highly  attenuative  materials[19]  by  analyzing  the  interference  pattern  between  the  in¬ 
cident  and  transmitted  continuous  ultrasonic  waves.  Our  technique  is  somewhat  similar  to  the  one 
proposed  for  the  measurement  of  attenuation  and  phase  velocity  in  liquids  by  Sedlacek  and  Asen- 
baum  [32].  They  obtained  an  interference  pattern  by  continuously  changing  the  distance  of  prop¬ 
agation  in  the  liquid  by  moving  one  of  the  transducers  either  closer  to  or  farther  from  the  other. 
Obviously,  this  can  be  done  for  propagation  in  a  liquid  but  not  in  a  solid. 

The  other  techniques  such  as  continuous  wave  (cw)  ultrasonic  spectroscopy  [  1 8],  for  simultaneous 
measurements  of  phase  velocity  and  attenuation  are  applicable  only  to  low-absoibing  materials. 
Standing  waves  are  set  up  in  the  plane  parallel  specimen  acting  as  a  resonant  cavity.  The  quality 
factor  Q  of  the  resonances,  defined  [33]  as  the  ratio  between  resonant  frequency  of  a  peak  and  its 
width,  gives  the  attenuation,  whereas  frequency  spacing  between  two  consecutive  resonances 
gives  the  phase  velocity.  This  technique,  however,  is  not  applicable  to  heterogeneous  materials  or 
to  materials  with  high  acoustic  attenuation.  In  fact,  in  these  cases,  multiple  resonances  are  either 
missing  or  are  very  broad  and  overlapping,  thus  rendering  the  measurements  meaningless,  partic¬ 
ularly  when  the  acoustic  parameters  must  be  carefully  evaluated  to  map  variations  in  material 
properties. 

These  above  difficulties  may  be  overcome  if  we  send  an  ultrasonic  cw  through  the  material  and 
cross-correlate  the  signal  from  the  transmitter  with  the  signal  from  the  receiver.  An  interference 
pattern,  which  contains  all  the  information  needed  for  the  evaluation  of  the  acoustical  parameters, 
can  be  obtained  by  continuously  varying  the  acoustical  path.  This  can  be  done,  in  the  simplest  pos¬ 
sible  way,  by  studying  a  material  specimen  cut  in  the  shape  of  a  wedge.  An  alternative  possibility, 
in  the  case  of  a  material  plate,  is  to  send  oblique  cw's  and  sweep  the  incidence  angle.  While  it  is 
obvious  that  plates  are  of  more  practical  interest  than  wedges,  we  prefer  to  devote  the  present  dis¬ 
cussion  to  the  analysis  of  wedges.  In  fact,  if  the  material  is  homogeneous,  the  only  difference  be¬ 
tween  waves  entering  wedges  at  different  locations  is  their  length  of  propagation  path  within  the 
material  (all  other  factors,  such  as  boundary  conditions,  being  identical);  therefore,  these  speci¬ 
mens  allow  us  to  concentrate  only  on  the  analysis  of  interference  effects,  without  having  to  worry 
about  the  angular  dependence  of  reflection,  refraction  and  mode  conversion  effects,  which  greatly 
complicate  the  physical  picture.  Consequently,  one  can  determine  absolute  attenuation,  since  the 
boundary  conditions  for  the  incoming  and  outgoing  cw's  are  the  same  at  any  given  point.  The  re¬ 
fraction  and  reflection  effects,  occurring  in  the  case  of  obliquely  incident  cw's  in  parallel  plates, 
will  be  considered  in  a  future  research  project. 

Besides  allowing  the  measurement  of  phase  velocity  and  attenuation,  our  technique  can  also  be 
used  to  detect  and  assess  quantitatively  material  flaws  and  inhomogeneities.  To  demonstrate  this, 
we  drilled  a  hole  in  one  of  the  specimens,  normal  to  the  direction  of  wave  propagation  and  scanned 
the  specimen  in  the  usual  fashion.  The  resulting  interference  pattern  clearly  shows  the  effect  of 
the  hole.  A  “signature  of  the  defect”  can  be  obtained  by  comparing  the  theoretical  acoustical  am- 
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plitude  curve,  expected  for  an  unblemished  specimen,  with  the  experimental  Our  results  sug¬ 

gest  the  possibility  of  an  approach  to  the  solution  of  the  inverse  problem  of  predicting  the 
characteristics  of  the  defect  through  an  analysis  of  the  experimentally  measured  acoustical  ampli¬ 
tudes. 

The  theory  and  experimental  setup  used  for  the  development  of  our  technique  are  explained  in 
3.2.1  and  3.2.2.  In  Section  3.2.3,  we  discuss  several  methods  for  the  analysis  of  the  experimental 
data  and  illustrate  them  with  results  obtained  for  an  unblemished  epoxy  specimen  and  for  the  same 
specimen  after  having  drilled  a  hole  in  it. 

3.2.1.  THEORY 

In  order  to  characterize  a  homogeneous  material,  we  fabricate  a  specimen  of  it  in  the  shape  of  a 
wedge  (Fig.  6)  and  place  it  in  an  ultrasonic  immersion  tank.  The  size  of  the  specimen  must  be  suf¬ 
ficiently  small,  so  that  it  slides  freely  between  the  transducers,  but  long  enough  so  that  a  large 
number  of  measurements  can  be  performed.  In  our  case,  the  distance  /  between  the  transducers 
was  about  3  cm  and  the  length  of  the  specimen  (x  range)  about  10  cm.  Continuous  longitudinal 
waves  are  propagated  through  the  sample  We  assume  that  the  wedge  angle  6  is  small  (in  our  case, 
e  =  13.8)  and  ignore,  for  the  time  being,  all  refraction  and  multiple  reflection  effects.  The  ultra¬ 
sonic  signals,  ST  sent  by  the  transmitter  T,  and  S„  received  by  the  transducer  R,  are  given  by  [12] 

ST  ~  FjA0cos  (<d/  +  4q)  ,  (3.2.1) 

SK  *  (l-y)  -o^yl 

+♦„),  (3.2.2) 

where  Ft  and  F„  are  the  response  factors  of  the  transmitter  and  of  the  receiver  and  associated  elec¬ 
tronics;  rl2  and  r2,  are  the  transmission  factors  between  medium  1  (water)  and  medium  2  (mate¬ 
rial)  and  vice  versa  (they  are  constant  since  both  the  wedge  angle  6  and  the  incidence  angle  are 
constant);  a, ,  a,,  v, ,  and  v2  are  the  attenuation  coefficients  and  phase  velocities  in  the  two  media; 
A0  is  the  amplitude  of  the  incident  cw;  /  is  the  (fixed)  distance  between  the  transducers;  and  y  is 
the  length  of  the  propagation  path  of  the  cw  in  the  material.  By  moving  the  wedge  (or  the  trans¬ 
ducers)  parallel  to  the  x  axis,  y  can  be  varied.  By  choosing  the  origin  of  the  x  axis  at  the  vertex  of 
the  wedge,  we  have 

ysx-tanO,  (3.2.3) 

where  8  is  the  angle  of  the  wedge. 

The  two  signals,  ST  and  S„ ,  are  correlated  through  an  electronic  mixer  (phase  detector).  The  cor¬ 
responding  non-normalized  cross-correlation  factor  is  given  by 
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(3.2.4) 


¥  «  JTrsr('>s* 

=Fj«Fji7'|27"jiAJmpI— ctj  (/-*y)  —  c^y]  3 , 

where 

3  ■  j^rco*(mf  +  4)  cos|o>0-^!-^)  +40J</r 


”  U-rcos[“(TT  +  ^)]‘" 


if'Tm  2“'*2»o-“(^+ * 


(3.2.5) 


The  cross-correlation  time  interval  2T  can  be  chosen  arbitrarily  large.  If  we  choose 


2T  »  — ,  (3.2.6) 

CD 

the  contributions  of  the  second  integral  in  Eq.  (3.2.5)  become  negligible  and 

3-rcos{a»[-^+(^)]}.  (3.2.7) 

It  then  follows  that 

y  =  y0  +  Cexp  (-ay)  cos  (by  +  4) ,  (3.2.8) 

where 

C  =  FTFRTnTixA\Texp  (-a,/) ,  (3.2.9a) 

a  =  Oj-a,,  (3.2.9b) 

b  »  co(j2-j,),  (3.2.9c) 

4  =  a )/s,,  (3.2.9d) 

and  j,  is  the  slowness  given  by 

j,  =  l/v  (i  =  1,2).  (3.2.10) 


In  Eq.  (3.2.8),  a  constant  term  y0  has  been  added  to  account  for  the  arbitrary  origin  of  the  exper¬ 
imental  acoustic  amplitude  scale. 

Upon  fitting  Eq.  (3.2.8)  to  the  experimental  data  (see  Sec.  3.2.3),  we  obtain  the  parameters  y0,  C, 
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a,  b,  and  ♦,  from  which  the  phase  velocity  v2  and  attenuation  a,  in  the  material  can  be  easily  de¬ 
duced,  since  the  corresponding  quantities  in  water,  v,  and  at,  are  well  known.  In  fact,  at  room 
temperature  293.5  K,  v,  =  1483  m/s,  and  since  a,  <<02,  a,  can  be  neglected.  It  follows  that 

otj  =  Qj+a-a, 


•wl  +  i)"-  (3.2.11.) 

Let  us  now  analyze  the  effect  of  refraction  (Fig.  7).  One  side  of  the  wedge  is  normal  to  the  incident 
wave  and,  therefore,  does  not  refract.  The  other  side  causes  a  deviation  0-6,  from  the  original 
propagation  direction,  where  0,  £  0  is  the  refraction  angle  at  the  wedge  surface  in  water  (we  are 
assuming  that  v2  >  v, ).  Consequently,  the  receiving  transducer  must  be  rotated  through  an  angle 
0  -  0J,  and  the  acoustic  path  increases  by  a  distance 


4  [isohg-'j 

-i(#-/0-y)  (O-©,)2-  (3.2.11b) 

However,  Eq.  (3.2.8)  still  holds,  but  the  values  of  the  parameters  C,  a,  b,  and  +  of  Eq.  (3.2.9) 
change  slightly.  In  particular,  when  refraction  effects  are  included, 

04  «  a+o, [l  +  (O-e*)2/^]  -a. 


b  [1  +  (0-e,)2/2]  -• 
v2  =  <  =  + - - - > 


(3.2.12) 


resulting  in  a  small  deviation  in  the  values  of  a2  and  v2.  Here,  0*  can  be  experimentally  measured 
or,  for  better  accuracy,  obtained  through  Snell's  law 

sinO^/v,  =  sin0/v2 

by  implementing  a  Newton  iteration  procedure  with  Eq.  (3.2.12). 

Finally,  we  discuss  briefly  the  effect  of  multiple  reflections.  Unless  r,2  =  r2,  =  1 ,  there  will  be 
multiple  reflections  of  the  waves  at  both  inside  surfaces  of  the  wedge  and  at  T  and  R.  A  number 
of  waves  are  thus  generated,  which  add  up  to  the  basic  signal  SK,  given  by  Eq.  (3.2.2).  Assuming 
e«l  and  considering  only  multiple  reflections  within  the  material  wedge,  we  note  that  if  the 
wave  is  reflected  n  times  by  each  inside  surface  of  the  wedge,  its  acoustic  path  within  the  wedge 
is  approximately  (2n+l)y.  Its  amplitude  is  multiplied  by  a  factor  R2',  where  R  is  the  reflection  co¬ 
efficient  of  the  material,  assumed  to  be  the  same  at  the  two  surfaces  of  the  wedge.  As  a  result,  we 
obtain  waves  with  cross-correlation  factors  still  given  by  Eq.  (3.2.8)  but  with  different  values  of 
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the  parameters  and  with  a  (slightly)  variable  transmission  coefficient  T2l .  The  interference  of  all 
these  waves  creates  a  slight  modulation  pattern,  which  needs  to  be  explicitly  taken  into  account 
for  a  better  fitting  of  the  experimental  curves.  Usually,  the  factor  R2"  and  the  high  attenuation  in 
the  material  reduce  the  amplitude  of  the  echoes  so  much  that  only  doubly  reflected  waves  (n=l) 
need  to  be  considered.  Also,  since  the  attenuation  increases  with  the  frequency,  the  modulation 
due  to  multiple  reflections  becomes  negligible  at  higher  frequencies. 

3.23.  EXPERIMENTAL  SETUP 

Before  starting  with  the  acquisition  of  the  experimental  data,  the  orientation  of  the  specimen  is 
adjusted  such  that  the  incident  beam  parallel  to  the  y  axis  scans  along  the  x  direction  in  the  wedge 
surface  parallel  to  the  x-y  plane.  Broad-band  transmitting  (T)  and  receiving  (R)  transducers  are 
mounted  on  the  same  bridge  and  scan  together  in  the  x  direction  and  index  in  the  z  direction.  Their 
spatial  positions  can  be  adjusted  independently.  They  are  aligned  using  pulse-echo  as  well  as 
pitch-catch  techniques  (Fig.  2).  The  position  and  orientation  of  the  transducer  R  are  adjusted  so 
that  the  received  signal  is  maximum. 

Once  the  orientations  of  specimen  and  transducers  are  fixed,  the  experimental  setup  is  switched 
to  the  configuration,  as  shown  in  Fig.  8.  Continuous  waves  of  suitable  frequency  (in  the  range 
2-20  MHz)  are  split  into  two  parts.  One  portion  of  these  waves  is  amplified  and  excites  a  broad¬ 
band  immersion  transducer  (T).  The  other  portion  is  fed  into  a  phase-sensitive  detector,  whose  dc 
output  (proportional  to  the  phase  difference  between  input  and  output  signals)  is  filtered,  digitized, 
and  transmitted  to  a  VAX  computer,  under  the  control  of  a  Tektronix  desk  top  computer.  The  fil¬ 
tered  signal  can  also  be  recorded  on  an  x-y  chart  recorder  as  a  function  of  the  propagation  path 
y  =  x  tan  9  in  the  material. 


3.23.  DATA  ANALYSIS  AND  DISCUSSION 

Typical  experimental  cross-correlation  curves  for  an  unblemished  epoxy  specimen  at  different 
frequencies  are  shown  in  Fig.  9.  They  clearly  show  the  exponential  decay  of  the  amplitude  and  the 
sinusoidal  oscillations,  as  predicted  by  Eq.  (3.2.8).  By  comparing  the  three  curves  obtained  at  dif¬ 
ferent  frequencies,  we  see  that  the  attenuation  o  increases  with  the  frequency,  as  expected.  We 
also  notice  that  the  intervals  between  successive  peaks  become  progressively  smaller  at  higher  fre¬ 
quencies,  as  predicted  by  Eq.  (3.2.9c),  which  shows  that  the  coefficient  b  of  Eq.  (3.2.8)  is  propor¬ 
tional  to  the  frequency. 

For  a  quantitative  fitting  of  the  experimental  data  and  consequent  extraction  of  the  parameters  of 
interest,  we  used  three  different  methods.  With  the  first  method,  we  obtain  the  five  parameters  y0, 
C,  a,  b,  and  $  by  fitting  five  points,  which  can  be  conveniently  chosen,  such  as  y,  =  0  (theoreti¬ 
cally,  at  the  edge  of  the  wedge  but,  in  practice,  sufficiently  apart  to  avoid  diffraction  from  the 
edge),  the  first  two  maxima,  first  minimum,  and  one  more  maximum  near  the  end  of  the  y  range 
(see  Fig.  10).  By  substituting  in  Eq.  (3.2.8)  and  noting  that  cos  {by  +  ♦)  is  equal  to  1  at  the  maxim? 
and  to  -1  at  the  minima,  one  easily  obtains 
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2k 


(3.2.13a) 


b  = 


(* -%>’ 


V0-  (¥2  +  2v3  +  ¥4)/4. 


o  *  tog [(v,- /(¥,-¥«,) ]/(>,- y2). 


(3.2.13b) 

(3.2.13c) 


+  »  atan  (- a/b )  -  by3  +  bjc,  n  =  0,  ±1 ,  ±2..., 


(3.2.13d) 


(¥!-¥«,) 

cos# 


(3.2. 13e) 


In  Eq.  (3.2.13d),  n  must  be  chosen  even  or  odd,  as  necessary  to  make  C,  in  Eq.  (3.2. 13e),  positive. 

Equations  (3.2.13)  and  (3.2.11a)  or  (3.2.11b)  allow  us  to  evaluate,  in  a  very  simple  fashion,  die 
attenuation  a,  and  phase  velocity  v2  in  the  material.  The  accuracy  may  not  be  too  good,  however, 
since  the  experimental  curve  is  determined  by  a  large  but  finite  number  of  points,  due  to  the  finite 
step  of  the  increment  in  the  x  direction  (0.025  mm).  Although  the  experimental  errors  are  rather 
small  (on  the  average,  much  less  than  1%),  it  may  be  difficult  to  pinpoint  the  precise  location  of 
the  maxima  (or  minima). 

A  second  method  for  the  analysis  of  the  experimental  data  consists  of  considering  all  the  maxima 
and  minima.  In  this  way,  it  is  possible  to  reduce  the  error  in  the  evaluation  of  the  interval  between 
two  successive  maxima  by  averaging  through  many  such  intervals.  This  method,  however,  re¬ 
quires  more  manipulation  of  data  and,  therefore,  is  much  more  time-consuming.  To  have  it  com¬ 
pletely  computerized  is  possible,  but  not  trivial,  since  local  random  fluctuations  may  be 
interpreted  as  maxima  or  minima,  thus  compromising  the  reliability  of  the  results.  This  method 
has  been  applied  for  the  data  analysis  contained  in  Ref.  19,  yielding  standard  deviation  for  the 
phase  velocity  of  the  order  of  0.5%. 

For  the  work,  we  adopted  a  different,  but  much  more  powerful  method,  which  consists  of  least- 
squares  fitting  all  the  (about  500)  experimental  points  to  the  theoretical  curve  given  by  Eq.  (3.2.8). 
In  this  way,  all  the  experimental  points  (rather  than  only  a  few  selected  ones)  contribute  to  the 
evaluation  of  the  parameters  of  interest.  A  computer  program  written  for  the  purpose  takes  only  a 
few  seconds  to  run  on  a  Digital  VAX  1 1/780  for  each  case  considered.  In  the  program,  a  dichot¬ 
omic  procedure  can  be  used  very  efficiently,  although  some  care  must  be  exercised  to  insure  con¬ 
vergence  to  the  global  minimum  (of  the  sum  of  the  deviations  square),  rather  than  to  a  local  one. 
With  this  method,  we  have  obtained  an  accuracy,  in  the  determination  of  v3,  of  better  than  0.1%. 
This  accuracy  could  be  further  increased  by  digitizing  the  experimental  acoustic  amplitude  mea¬ 
surements  over  a  wider  range  (in  our  case,  the  range  was  0-255),  thus  reducing  accordingly  the 
experimental  errors. 

Although,  in  principle,  the  same  accuracy  can  also  be  achieved  in  the  evaluation  of  the  attenua¬ 
tion,  in  practice  the  standard  deviations  here  are  much  larger,  because  of  possible  losses  due  to 
surface  roughness,  cutting  of  the  wedge  surface  that  is  not  perfectly  straight,  and  local  inhomoge- 
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neities.  On  the  other  hand,  the  attenuation  measurement  can  be  used  to  detect  these  or  other  irreg¬ 
ularities.  It  should  also  be  pointed  out  that,  unless  they  are  properly  taken  into  account,  (multiple) 
internal  reflections  can  be  a  large  source  of  error  in  the  determination  of  the  attenuation,  especially 
at  low  frequencies. 

In  order  to  demonstrate  the  applicability  of  our  technique  to  the  detection  of  defects,  we  have 
drilled  a  hole  (of  0.25-cm  diam.)  in  one  of  our  epoxy  specimens,  normal  to  the  cw  propagation 
direction,  and  measured  the  corresponding  acoustical  amplitude  pattern.  The  result  is  shown  in 
Fig.  1 1 .  By  comparing  Fig.  1 1  with  Fig.  9(a),  we  can  very  clearly  observe  the  effect  of  the  hole  in 
the  region  between  y  =  0.3  and  0.5  cm.  The  effect  can  be  best  observed  by  subtracting  the  exper¬ 
imental  acoustic  amplitude  from  the  theoretical  least-squares  fitted  curve.  For  an  even  better  rep¬ 
resentation,  after  having  determined  where  the  defect  lies,  one  can  repeat  the  least-squares  fitting 
excluding  the  range  affected  by  the  defect.  In  this  way,  one  obtains  values  of  04  and  v2  that  are 
not  affected  by  the  presence  of  the  defect.  Also,  by  subtracting  the  experimental  data,  one  obtains 
the  true  “signature”  of  the  defect  [Fig.  12(b)].  For  comparison.  Fig.  12(a)  shows  the  random  fluc¬ 
tuations  obtained  when  the  defect  is  absent.  Other  smaller  holes  (down  to  0.39  mm)  were  drilled 
elsewhere  in  the  specimen;  their  signatures  were  noticeable,  although  much  less  conspicuous.  If 
better  precision  of  measurements  were  available,  we  could  have  better  representations  of  these  de¬ 
fects  and  consequently  reduce  the  background  noise.  By  accumulating  results  like  Fig.  12(b)  for 
different  kinds  and  sizes  of  defects,  one  can  learn  to  solve  the  inverse  problem  of  detecting  and 
quantitatively  assessing  an  unknown  defect. 

Beside  the  information  contained  in  the  envelope  of  the  curve  of  Fig.  12(b),  the  periodicity  of  the 
oscillations  and  the  phase  shift  found  through  the  defect  range  can  provide  information  about 
phase  velocity  changes  in  the  defect  region.  In  our  specimen,  we  have  found  a  decrease,  albeit 
very  small,  in  the  phase  velocity,  possibly  due  to  the  acoustoelastic  effect  [34]  of  the  residual 
stress  introduced  by  the  drilling  of  the  hole.  The  effect,  however,  was  too  small  (only  slightly  larg¬ 
er  than  the  standard  deviation  for  v2)  to  justify  a  quantitative  analysis. 


3 3  ULTRASONIC  SPECTRAL  ANALYSIS 

The  pulse-echo  [35]  technique  which  is  commonly  used  in  measurements  of  attenuation  is 
inadequate,  due  to  lack  of  multiple  echoes,  as  seen  in  acoustically  transparent  parallel  plate  metallic 
specimens.  In  the  previous  section  we  used  wedge-shaped  specimens  to  measure  attenuation  in 
these  materials,  using  an  ultrasonic  interferometric  technique.  Such  a  technique,  however,  requires 
a  special  geometry  of  the  specimen  and  therefore  can  be  of  limited  use.  The  attenuation  in  a 
material  is  loss  in  energy  due  to  absorption  and  scattering.  In  general,  as  pointed  out  earlier,  the 
signal-to-noise  ratio  is  very  poor  and,  depending  on  the  thickness  of  the  specimen,  one  may  not  be 
able  to  observe  a  transmitted  pulse.  One  can  determine  the  attenuation  from  this  pulse  using  a 
substitution  method  [36]  or  by  time  delay  spectroscopy  [37].  Neither  of  these  methods  yields 
information  about  the  frequency  dependence  of  the  attenuation,  unless  by  sweeping  the  frequency. 

The  attenuation  of  propagating  ultrasound  through  heterogeneous  materials  such  as  fiber- 
reinforced  epoxy  composites  containing  chopped  fibers  with  random  orientation  is  predominantly 
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due  to  scattering  from  the  fibers  as  well  as  from  discontinuities  such  as  pores  and  is  consequently 
frequency  dependent.  In  this  section,  we  determine  the  frequency  dependent  attenuation  from  the 
downshifted  spectra  of  the  transmitted  broadband  wavepacket.  Such  a  technique,  first  used  for 
medical  ultrasound,  can  be  very  useful  for  characterization  of  acoustically  very  lossy  composites, 
because  spectral  shifts  for  such  materials  can  be  particularly  significant.  It  also  allows  use  of 
parallel  plate  specimens.  Analytical  models  and  experimental  results  are  presented  in  sections 
(3.3.1)  and  (3.3.2). 

33.1  THEORY 

An  ultrasonic  wave  propagating  through  an  inhomogeneous  material  suffers  losses  in  energy  due 
to  scattering,  which  can  be  caused  by  spatial  variations  in  absorption  properties  or  in  material 
density  and/or  in  wave  velocity.  We  assume  that  the  attenuation  in  a  signal  propagating  through 
such  a  dispersive  material  is  frequency  dependent  and  is  given  by 

«(/)  =  Uff  (3.3.1) 

where  a,,  and  n  are  material  characteristics. 

We  briefly  review  and  extend  the  mathematical  analysis  of  Narayana  and  Ophir  [38].  Let  us 
consider  a  rectangular  burst  of  amplitude  A0,  duration  T0  and  center  frequency  /„,  incident  on  a 
plane  parallel  plate  of  composite  material.  Such  a  pulse  in  the  time  domain  can  be  written 

h  (r)  =  A0cos  (2nf0t)  for  111  <  r0/2  (3.3.2) 

=  0  for  Id  >  r0/2 
Its  Fourier  transform  [39]  is  given  by 

H[f)  =  A\T0  IQ  (f+f0)  +  G  (/-/„)] 


where  Q  (/) 


sin  (nTff) 
nTJ 


The  frequency  spectrum  of  the  incident  pulse  is  then  given  by 


(3.3.3) 


Agr0sinrcr0(/-/0)  (334) 

w  ‘  *r0  (/•-/„) 

When  such  a  pulse  is  transmitted  through  an  attenuating  material,  both  the  amplitude  of  its  center 
peak  and  its  bandwidth  decrease.  Further,  if  the  attenuation  is  frequency  dependent  (assume  n  >  1), 
then  the  amplitude  at  higher  frequencies  is  attenuated  more  than  at  lower  frequencies. 
Consequently  there  is  an  apparent  frequency  shift  of  the  peak.  The  attenuation  of  the  material  can 
be  determined  from  this  shift  in  frequency  of  the  central  maximum  (peak). 


13 


We  describe  the  medium  of  propagation  by 


(3.3.5) 


where  a  (/)  is  given  by  Eq.  (3.3.1)  and  z  is  path  of  propagation  in  the  lossy  material.  The  spectrum 
of  the  signal  transmitted  through  the  material  is  given  by 


T(f)  =I(J)®AUf) 


(3.3.6) 


In  order  to  find  the  shift  in  center  frequency  /0,  we  set 


dT(f) 

df  /«,. 


(3.3.7) 


where  fe  is  the  center  frequency  of  the  received  wavepacket,  i.e. 


d  r ((-/o) 

df[  2  *T0(f-f0)  *  \JmJt  = 


(3.3.8) 


Upon  simplification,  one  obtains. 


coskt0  ( fc  -/0)  -  S“rJcfl/or  Ina°^~ 1 V'  ~/o)  + 1 1  "  0 


(3.3.9) 


x  =  (4-/0)kT0 


(3.3.10) 


then  Eq.  (3.3.9)  can  be  written  as 


«v:-  -  — T-J—n 

”c  z  Ltan  x  xJ 


(3.3.11) 


r-i — ii 

Ltan  x  xJ 


(3.3.12) 
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(3-3.13) 


««o£ 


*  -YTn 


or 


K  YTa 


ou  * - -  nepers  cm  MHz 

^  znfe~l 


(3.3.14) 


For  the  purpose  of  analyzing  the  experimental  data,  one  can  write  Eq.  (3.3.14)  in  a  linear  form  by 
taking  the  logarithm  of  both  sides.  From  this  linear  relationship,  one  can  determine  n  or  o^,  by 
measuring  fe,  for  two  different  values  of  T0  or  /0.  The  attenuation  in  the  material  at  frequency  /0  is 
then  determined  by 

a  =  8.6860^3  dB  cm~l  (3.3.15) 

where  the  material  parameters  a0  and  n  have  been  determined  by  the  procedure  outlined  above. 
33.2  EXPERIMENTAL  SET-UP  AND  RESULTS 

Figure  13  shows  the  experimental  set  up.  Two  matched  2.25  MHz  broadband  immersion 
transducers  were  used  to  transmit  and  receive  a  rectangular  toneburst  through  a  parallel  plate 
specimen  (70  mm  x  40  mm  x  8  mm).  A  rectangular  toneburst  of  a  given  frequency  is  generated  by 
gating  a  continuous  wave  rf  signal  with  a  very  precise  rectangular  pulse.  The  amplified  toneburst 
is  used  to  excite  the  transmitter.  Both  the  incident  and  the  received  signals  are  monitored  on  a  scope 
to  make  sure  that  the  amplifier  gains  are  proper  and  the  shape  of  these  pulses  is  rectangular  as 
desired. 

First,  the  transducers  are  aligned  without  the  specimen  and  the  spectrum  of  the  pulse  propagating 
through  water  is  analyzed.  The  effect  of  amplifier  gains  on  the  shape  and  peak  frequency  of  the 
pulse  transmitted  through  water  is  studied  by  attenuating  the  incident  pulse  (Fig.  14).  It  is  clear  that 
the  peak  frequency  does  not  shift  as  a  function  of  artifacts  such  as  amplifier  gain  etc.  It  is  expected 
this  will  be  the  case  even  in  the  presence  of  the  specimen. 

Figures  15  (a-c)  show  the  changes  to  the  center  frequency  of  the  peak  for  various  durations  of  the 
incident  rectangular  toneburst.  These  effects  are  caused  by  attenuation  in  the  specimen.  For  each 
T0,  the  upper  trace  shows  the  spectrum  of  the  pulse  propagated  through  water.  When  the  specimen 
is  inserted  between  the  transducers,  three  things  happen:  (1)  the  amplitude  of  the  peak  drops,  (2) 
the  position  of  the  center  peak  is  downshifted,  and  (3)  the  bandwidth  of  the  center  peak  decreases. 
The  incident  pulse  for  water  as  specimen  is  attenuated  by  20  dB,  to  keep  its  spectrum  on  the  page 
along  with  that  of  the  composite  material  for  comparison  and  analysis.  As  r0  increases,  e.g.  from 
broadband  to  narrowband  (quasi-condnuous  waves),  the  shift  in  the  frequency  decreases,  as  one 
would  expect  qualitatively.  But  the  number  of  sideiobes  also  increases  as  T0  increases  from  6  usee 
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to  1 1  nsec.  Using  these  measured  values  of  f0  and  fc  in  Eq.  (3.3.14),  we  find  from  a  set  of 
simultaneous  equations,  a  value  of  n  =  58,  which  is  rather  high.  But  it  indicates  that,  due  to  various 
heterogeneities  present  in  the  specimen,  the  attenuation  is  mostly  due  to  scattering  and  depends 
very  strongly  upon  frequency.  The  amount  of  attenuation  due  to  scattering  at  2.7  MHz  was 
determined  to  be  about  16  to  20  dB/cm  depending  upon  the  location  in  the  specimen.  The  total 
attenuation  in  these  materials  determined  earlier  [40]  by  using  a  wedge  specimen  was  found  to  be 
~32  dB/cm.  From  the  above  data,  using  a  substitution  technique  analysis,  the  total  attenuation  turns 
out  to  be  27  dB/cm,  which  agrees  well  with  previous  results.  This  variation  of  15%  in  the  measured 
value  can  easily  be  caused  by  heterogeneities  in  the  material.  By  measuring  the  frequency 
dependent  component  of  the  attenuation,  it  is  demonstrated  that  in  these  materials  the  attenuation 
is  indeed  dominated  by  losses  in  energy  due  to  scattering. 

4.  SUMMARY 

In  this  report  we  have  shown  how  to  characterize  highly  attenuative  materials  which  can  barely 
transmit  ultrasound.  The  type  of  energy  used  for  nondestructive  evaluation  of  thick  composites 
must  be  able  to  penetrate  through  the  thickness  of  the  specimen  under  interrogation.  One  needs  a 
sensitive  receiver  and  good  signal  processing  for  extracting  signals  due  to  a  flaw.  A  specimen 
which  is  “thick”  (i.e.,  poor  signal/noise  ratio)  for  one  type  of  NDE  method  may  not  be  “thick”  for 
another  NDE  method,  provided  the  method  is  capable  of  evaluating  the  specific  flaws  encountered. 

We  have  demonstrated  that  it  is  possible  to  measure  the  attenuation  accurately  in  highly  absorbing 
materials  by  using  a  transmission  technique  and  a  properly  designed  specimen.  The  internal 
composition  of  the  material  is  clearly  depicted  by  the  attenuation  curve. 

We  have  discussed  the  application  of  ultrasonics  to  characterize  thick  composites.  Ultrasonic 
techniques  based  on  propagation  of  elastic  waves  through  the  material  are  generally  applicable  to 
all  types  of  highly  attenuative  materials  whereas  eddy  current  techniques  are  applicable  only  to 
electrically  conducting  composites  such  as  graphite-epoxy.  The  data  acquired  for  the  various  types 
of  flaws  can  be  transformed  into  computer-generated  images  of  spatial  variation  of  selected 
material  properly. 

We  proposed  an  ultrasonic  technique  for  the  simultaneous  determination  of  the  phase  velocity  and 
attenuation,  which  can  be  used  for  the  characterization  of  highly  attenuative  materials.  The  tech¬ 
nique  is  based  on  the  interference  between  the  incident  (reference)  cw  signal  and  the  signal  after 
propagation  through  the  material.  The  interference  pattern,  obtained  when  the  acoustical  path  is 
continuously  varied,  contains  all  the  information  necessary  to  evaluate  the  phase  velocity  and  at¬ 
tenuation  of  the  material  and  any  eventual  anomaly  or  irregularity  in  the  material.  In  order  to  ob¬ 
tain  a  continuous  variation  of  acoustical  path  without  the  complications  arising  from  changing 
boundary  conditions,  we  fabricated  a  material  specimen  in  the  shape  of  a  wedge  and  scanned  it  in 
the  direction  normal  to  its  edge.  Characteristic  interference  patterns  have  been  obtained  from  the 
acoustical  amplitude  at  different  frequencies.  Several  methods  were  adopted  for  the  analysis  of 
the  experimental  data.  We  have  found  that  a  least-squares  Fit  to  the  theoretically  predicted  curve 
can  be  a  very  efficient  procedure,  which  allows  the  determination  of  the  phase  velocity  with  an 
accuracy  of  better  than  0.1%.  Larger  standard  deviations  are  found  for  the  attenuation,  but  they 
reflect  physical  characteristics  of  the  specimen,  such  as  inhomogeneities  or  imperfect  surface.  By 
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drilling  a  hole  through  the  specimen,  normal  to  the  direction  of  propagation  of  the  cw,  we  simu¬ 
lated  the  presence  of  a  defect  and  showed  that  our  technique  is  quite  capable  of  detecting  it.  By 
subtracting  the  experimental  acoustic  amplitudes  from  the  theoretical  curve,  we  have  shown  that 
it  is  possible  to  obtain  a  distinctive  track  of  the  defect,  which  can  be  used  for  its  quantitative  as¬ 
sessment,  once  a  sufficient  number  of  signatures  of  defects  of  different  size  and  shape  have  been 
collected  as  reference. 

Measurement  of  attenuation  by  frequency  shift  in  the  spectra  of  a  transmitted  wavepack-t  of 
known  shape  allows  use  of  parallel  plate  specimens.  It  is  a  viable  technique  for  characterization 
of  such  materials.  Analysis  and  technique  could  be  extended  to  other  types  of  wavepacket  and 
specimens  of  other  composite  materials,  in  which  scattering  of  the  ultrasound  is  the  dominant 
cause  of  attenuation.  The  model  should  be  extended  to  correlate  the  frequency  shift  in  the  side 
lobes  of  a  transmitted  pulse.  It  is  expected  that  side  lobes  will  have  higher  frequency  shifts  and 
may  give  increased  accuracy. 
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FIG.  1.  Two  wedge  specimens  (epoxy  and  composite  material)  and  a  step 
specimen  (composite  material). 
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FIG.  3.  Attenuation  curve  for  the  step  specimen. 
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FIG.  4a.  Attenuation  curve  for  the  epoxy  wedge  specimen. 
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FIG.  4b.  Attenuation  curve  for  the  epoxy  wedge  specimen  with  a  flaw. 


TNoIumm  (cm) 

FIG.  5b.  Attenuation  curve  for  the  flawed  regions  in  the  composite  wedge  specimen. 
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FIG.  5a.  Attenuation  curve  for  the  unflawed  regions  in  the  composite  wedge  specimen 
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FIG.  6.  Arrangement  of  specimen  and  transducers  (refraction  effects  ignored). 
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FIG.  7.  Arrangement  of  specimen  and  transducers  (refraction  effects  included). 
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FIG.  8.  Experimental  setup  for  ultrasonic  interferometric  technique. 
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FIG.  9.  Typical  experimental  cross-correlation  curves  for  a  specimen  without  defect  at  different 
frequencies:  (a)  1 1  MHz,  (b)  15  MHz,  and  (c)  19  MHz. 
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FIG.  10.  Graphical  representation  of  the  five-point  data  analysis. 


31 


ATTENUATION 


FIG.  11.  Cross-correlation  curve  for  a  specimen  with  a  hole  drilled  in  it. 
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FIG.  12.  Difference  between  the  theoretical  least-squares  fitted  curve  and  the  corresponding 

experi  mental  curve  (a)  for  a  specimen  without  defect  and  (b)  for  a  specimen  with  defect. 
The  signature  of  the  defect  is  cleaHy  visible. 
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FIG.  14.  spectrum  of  toneburst  propagated  through  water. 
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FIG.  15a.  Spectrum  of  toneburst  of  duration  6  ps  propagated  through  composite  and  water. 
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FIG.  15b.  Spectrum  of  toneburst  of  duration  7.2  ps  propagated  through  composite  and  water. 
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FIG.  15c.  Spectrum  of  toneburst  of  width  lips  propagated  through  composite  and  water 
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